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Abstract 

Context. We present a new search for variable stars in the Galactic globular cluster M28 (NGC 6626). 

Aims. The search is based on a series of BVI images obtained with the SMARTS Consortium's 1.3m telescope at Cerro Tololo Inter- 
American Observatory, Chile. 

Methods. The search was carried out using the ISIS v2.2 image subtraction package. 

Results. We find a total of 25 variable stars in the field of the cluster, 9 being new discoveries. Of the newly found variables, 1 is an 
ab-type RR Lyrae star, 6 are c-type RR Lyrae, and 2 are long-period/semi-regular variables. V22, previously classified as a type II 
Cepheid, appears as a bona-fide RRc in our data. In turn, V20, previously classified as an ab-type RR Lyrae, could not be properly 
phased with any reasonable period. 

Conclusions. The properties of the ab-type RR Lyrae stars in M28 appear most consistent with an Oosterhoff'-intermediate classi- 
fication, which is unusual for bona-fide Galactic globulars clusters. However, the cluster's c-type variables do not clearly support 
such an Oosterhoff" type, and a hybrid Oosterhoff" I/II system is accordingly another possibility, thus raising the intriguing possibility 
of multiple populations being present in M28. Coordinates, periods, and light curves in diff'erential fluxes are provided for all the 
detected variables. 
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1. Introduction 

M28 (NGC 6266) is a moderately reddened [E(B - V) = 
0.40 mag] and bright (My = -8.16 mag) globular cluster (GC) 
located at a low Galactic latitude (b = -5. 58)[^ Though rela- 
tively close, at a distance from the Sun of only 5.5 kpc, it re- 
mains a relatively ill- studied cluster, likely due to the unfavor- 
able (and highly variable) foreground reddening. Indeed, vari- 
ability studies of the cluster have so far been mostly restricted 
to photographic data ( ^Wehlau & Butterworthnl99Q i Wehlau & 
|Clement|1990| ), and modem color-magnitude dia grams (CMDs) 
have only recently bee n presen ted in the literature ( Davidge et al. 
1996t [Rosenberg et al. 2000; Testa etal. 2 001 ; Alonso-Garcia et 



al.|2012| ). 



Still, M28 appears as a particularly interesting object for at 
least three reasons: first, it has been found to have a disk-like 
orbit ( Rees & C udworth|199H[Cudworth & Hanson] 1993) , thus 
making it one of the most metal-poo r members of the so-called 
"thick disk" family of GCs (see, e.g.,"Armandrofl r & Zinn|1988l ). 
Second, among clusters of similar metallicity ([Ee/H] ^ -1.32), 
M28 stands out as having a horizontal branch ( HB) morph ology 
strongly skewed towards the blue (see also | Alcaino| 198 1 1 ), thus 



* Based on observations obtained with the SMARTS Consortium 
1.3m telescope at the Cerro Tololo Inter- American Observatory, Chile 

^ Unless ot herwise noted, all cluster parameters in this paper are from 
Harris| ( [T996l Dec. 2010 update). 



making it a "second parameter" cluster. The latter, if interpreted 
in terms of age, would accordingly indicate that the metal-poor 
tail of the thick disk is at least as old as the oldest components of 
the Galactic halo. Third, the cluster has been classified into an 
QosterhoflFl ( |1939| |1944 ) I (Ool) type, which is consistent with 
the relatively high metallicity of the cluster but in conflict with 
its HB type: with rare exceptions, blue HB clusters are com- 
monly associated wit h type Poll ([Contreras et al.|20Q5||201Q| ). 



According to the 



Clement et al. 



( [2001t cataiogpl there are 



at present 24 variable stars known in the field of the cluster (in 
addition to a millisecond pulsar; Ly ne et al.|1987[ ). Of these, VI- 
VI 6 were discovered by Sawyer (19 49[), whereas V17-V2 4 were 



first discovered/reported on by , Wehlau & Sawyer Hogg ( 1982 



. Updated ephemerides for t hese stars were provided in 
Wehlau & Butterworth[ ([T99Q') and " Wehlau & Clement[ ( [T990| ), 
which remain the most recent papers to deal with the (photo- 
graphic) light curv es of M28 variable stars in a systematic way. 
Rees & Cudworth ([1991) reported, again based on photographic 
plates, on three additional candidate variables in the direction 
of the cluster, which however have not yet been systematically 
studied, or eve n incorporated into the electronic version of the 
Clement et al.[ catalog. 

There are strong reasons to believe that the current variable 
star tally for M28 is incomplete. Eirst, while M28 is a very 



http : //www . astro . utoronto . ca/~cclement/read . html 
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concentrated cluster, the photographic material which has been 
used in all previous variability studies does not allow the cluster 
core to be reliably res olved. S econd, mo dem image- subtraction 
techniques (e.g., Ala rd & Lup ton 1998J |Alard|[IDQQt |Bramich| 
|20Q8| ), when applied to modem CCD images, have recently 
been providing rich harvests of variable stars towards the cen- 
ters of even the previously best stu died globular clusters (see, 
e.g., Contreras et al.|2QlQ Kuehn et al. 2011 , for recent exam- 
ples). Accordingly, the main purpose of this paper is to provide 
the first CCD-based time-series photometry for the central re- 
gions of M28, in order to search for additional variables that may 
have gone unnoticed in previous studies and to provided updated 
periods and light curves for previously studied variables. 

We begin in §2 by describing out dataset and reduction tech- 
niques. In §3, we describe our variability results. A summary of 
our results is finally provided in §4. 



2. Observations and data reduction 

Our variability search is based on a total of 39 images in each of 
the Johnson-Cousins BVI filters, obtained in service mode with 
the SMARTS 1.3m telescope at CTIO. The exposure times were 
70, 25, and 15 seconds in B, V, and /, respectively. The images 
cover the period between June 2 and July 24, 2004, with 6 to 10 
images obtained every 3 to 5 days. 

The images were reduced in the standard way, using IRAF 
tasks p] A search for stel lar variability was performed using ISIS 
v2.2 ( |Alard & Lupton|[l998t [Alard|[20Q0| ), but modifying the 
search module so as to include candidates from all the avail- 
able filters. For the variable star candidates detected, periods of 
variabi lity were obtained using phase dispersion minimization 
(PDM; [StelTingwerfl 1978) . 



3. Variable stars 

Based on the variability search conducted with ISIS, we were 
able to find 25 variable stars with good photometry to reveal 
clean light curves. A finding chart with the new discoveries is 
given in Figure [T] (For clarity, previously known variables are 
not shown in this chart. Finding charts indicating the positions of 
those variables, including those that fall outside our field of view, 
are provided by We hlau & Sawyer Ho gg 1984 ) Light curves 
for the individual variable stars (in relative flux units) are pro- 
vided in Figures [2] to |5] (^/^c^r(9^/c version only). The detected 
variables include 6 ab-type RR Lyrae, only one of which is a 
new discovery. In addition, 9 c-type RR Lyrae were found, 6 of 
which being new discoveries, of which one is a possible field in- 
terloper. The previously reported RR Lyrae stars are confirmed 
in our data, as are the 4 type II Cepheids. However, two of these 
candidate Cepheids could not be cleanly classified on the basis 
of our data. In particular, for V21 we are not able to phase the 
data properly, which might suggest tha t the star is a long-period 
variable (LPV). This is surprising, since ] Wehlau & Sawyer Hogg 
(1984) present a nicely phased light curve for this star, using a 
per iod of 29.93 d. Unfor tunately, the star appears saturated in 
2012) study; however, its position in 



the Alonso-Garcia et al. 



our own preliminary CMD (not shown) is not inconsistent with 
an LPV classification. For V22, instead of the favored period of 



0.99538 d in the 2009 versio n of the [Clement et al.|(|200T]) cat - 
alog, which also comes from 'Wehlau & Sawyer Hogg| ( [1984 ), 
our data favor a much shorter period, of 0.323 d, and so this 
star may be better classified as a c-type RR Lyrae. The position 
of the star in the Alon so-Garcia et al.jCMDs is fully con sistent 
with a c-type classification. In [Wehlau & Sawyer Hogg[ ( [1982| ) 
that reported the discovery o f V22, a period of 0.498828 d was 
originally given, whereas in Wehlau & Sawyer Hogg (1984) a 
period of 0.498 d was given as an alternative, with a note also 
emphasizing that the star is blended in their images, which is 
supported by the large amount of scatter in their derived light 
curve. V21 and V22 are both fairly close to the center for such 
early photographic work. 

The case of V20 is particularly intriguing. While previously 
a period of 0.49774 d was reported ( [Wehlau & Sawyer Hogg| 
1984 ), and the star classified as an ab-type RR Lyrae, we could 
not phase our data properly with this period. A period around 
0.4948 d provides a better fit to our data, though with signif- 



icant scatter. While this is consistent with [Wehlau & Sawyer 
"Hogg*s comment that the star is blended, its mean magnitude in 
our preliminary photometry is not noticeably brighter than the 
magnitudes of other RR Lyrae stars in the cluster, and neither is 
it much redder than many of the other RR Lyrae. The light curve 
shape supports an ab-type classification, even though a cleaner 
light curve for a definitive assessment would be desirable. The 
star's position in the Al onso-Garcia et al.[ (|2012|) CMDs is also 
fully consistent with an RRab classification. Unfortunately, our 
data set is insufficient for an assessment of the incidence of the 
Blaz hko[ ( |l907 ) eff'ect among our studied stars. 



As to the remaining stars, namely V2, V3, V7, and VIO, we 
confirm the results of previous studies, in that the stars appear 
to vary on long timescales. For V2 and V3, we provide periods 
for the first time. V7 had a previously reported period of 320 d, 
which is much too long a timescale for our dataset to provide 
meaningful in formation regarding p eriodicity. In any case, its 
position in the [Alonso-Garcia et aT] ( ^2012) CMDs clearly con- 
firm that it is a fairly metal-rich red giant star in the background. 
Finally, for VIO we are again unable to phase the light curve 
with a significant period. The remaining known variables not 
in our study generally fall outside our field of view: VI, V8, 
V9 (an ab-type RR Lyrae), V6 (an LPV), V15, V16 (both ab- 
type RR Lyrae belonging to the field). V 14 is likely not variable, 



according to the Clement et al. 


([2001 ) catalog (see Wehlau &[ 


Butterworth 1990 ), even though 


Wehlau & Sawyer Hogg (1984j) 



provide possible periods around 0.2688 d for this star. 

Most of the variable star candidates listed by [Rees & 
Cud worth] ( 1991| ) could not be confirmed in our study. However, 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



based on the position in their finding chart, and as indicated in 
Table [3] RC133 may be the same star as our NV5, even though 
the coordinates published in their paper do not exactly match 
ours. The matching to (Rees & Cudworth| was done by triangula- 
tion and by eye. 

Table [3] summarizes our results. Column 1 gives the star 
name, whereas columns 2 and 3 give their right ascension and 
declination, respectively (epoch J2000). We derived celestial co- 
ordinates for all stars detected in our photometric study by com- 
parison with bright stars obtained from the Two Micron All 
Sky Survey (2MASS; Skrut sE^et al.|2006| ) catalog stars avail- 
able through the Infrared Processing and Analysis Center (IR\C) 
website. Around 1300 stars were available as astrometric refer- 
ences within our image. A third-order polynomial fit, done with 
the IRAF task mscpeak, produced dispersions of cr ~ 0.2'', con- 
sistent with the catalog precision. Column 4 gives our best-fitting 
period, with as many decimals places as justified on the basis of 
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Figure 1. Finding chart for the newly discovered variable stars in M28. The large image corresponds to a V-band image of the 
cluster, with the scale of the drawn blue square being 6.5 x 6.5 arcmin^. North is up and East is to the left. The five inner regions of 
the cluster where new variables have been detected are shown as blue circles superimposed on the main image. Blow-ups are shown 
for each of these regions, based on an F555W (V-band) Hubble Space Telescope image obtained with the WFPC2 camera (ID 6625, 
PI R. Buonanno). Each such circle covers a region of 30 arcsec diameter. In the top right, a 7.5 x 7.5 arcsec^ subregion of the circle 
containing NV3, NV6, and NV7 is shown on an expanded scale, for clarity. 



our data only, whereas column 5 gives the previously favored 
period. Column 6 gives the variability type favored by our data, 
whereas column 7 gives the previously favored variability type. 
The final column gives an alternative identification for the star, if 
available. Light curves for the individual variable stars (in rela- 
tive flux units) are provided in Figures [2] to |5](^/^c^r(9^/c version 
only). 

Among the ab-type RR Lyrae in M28, based on the magni- 
tude values published in the online [Clement et aL] ( [2QQl ) cat- 
alog, we associate as likely field variables not only stars VI 5, 
VI 6, and V24 (all of which are background stars), but also V9, 
which is about 1 mag brighter than the other bona-fide cluster 
RR Lyrae stars, and is thus likely a foreground RR Lyrae. The 
newly discovered RRab, NVl, appears to be a bona-fide cluster 
member, judging from its position in the Alonso-Garcia et al. 
( |2012j ) CMDs. 



As far as c-type RR Lyrae are concerned, no previous field 
interlopers had been known. Among the newly detected RRc's, 
only NV3 appears clearly to be a foreground RR Lyrae, being as 
it is about 2 mag brigh ter than other RR Lyrae stars in the cluster, 
again according to the [Alonso-Garcia et al.| ( |2Q12| CMDs. 

For the RR Lyrae stars that are clearly cluster members, one 
obtains the following average quantities and population ratios: 
<Pab) = 0.602 d (or 0.614 d, if we drop V20); <Pc) = 0.304 d; 

= NJNc+ah = 0.44 (or 0.47, if V20 is excluded). In addition, 
the minimum ab-type RR Lyrae period (/^ab,min) and maximum 
c-type RR Lyrae period (/^c,max) are 0.4916 d and 0.329 d, re- 
spectively. These values are quite unusual among known GCs. In 
particular, both (Pah) and /^a b,min are more typical of Qosterhoff - 
intermediate systems (e.g., Catelan""2009[ Catelan et al.[[2012 ). 
As discussed by Catelan et al. ( 2012 ), these are the two quantities 
most strongly defining Oosterhoff' status. Still, M28's (Pc) value 
is clearly more typical of Oosterhoff' type I systems, whereas 
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the high c-type number fraction is more common in OosterhofF 
type II systems. A "hybrid" OosterhofF I/II system is accord- 
ingly another possibiHty, which raises the question of whether 
there might be multiple populations in this fairly massive clus- 
ter. Given the known connections between Oosterhoff status and 
formation history (e.g.,|Catela n 2009 ; C lementini|2010[ van den 



Bergh'2011! Catelan et al. 2012 ), it would certainly be of interest 
to study the variable star content in M28 using a more extensive 
dataset. 

To close , we comment on the status of NV8. The star's po- 
sition in the I Alonso-Garcia et al. ( 2012| ) CMDs, fairly close to 
the tip of the red giant branch and at a fairly red color, is consis- 
tent with an LPV classification. This is why this star is grouped 
with other LPVs in Figure [5] On the other hand, according to 
our time-series data, a period around 13.53 d provides a decent 
match of the light curve (Fig. [5]). Such a period is, however, more 
consistent with a W Virginis (type II Cepheid) classification. 
Further analysis will clearly be necessary, before the variability 
nature of this star can be conclusively established. 
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4. Summary 

We have presented the results of a new search for variable stars 
in M28. Our search has led to the discovery of a number of pre- 
viously unknown variables, most of which are c-type RR Lyrae 
stars. The properties of the ab-type RR Lyrae stars are most con- 
sistent with an Oosterhoff-intermediate classification, but this is 
not clearly supported by the properties of the c-type RR Lyrae 
stars. A "hybrid" Oosterhoff' I/II classification is thus possible, 
raising the question as to whether multiple populations may be 
present in this fairly massive cluster. More extensive, higher- 
quality datasets will be required to put the properties of this clus- 
ter on a firmer basis. 
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Table 1. Variable stars in M28 



Name 


RA (J2000) 


DEC 


Period 


Old period 


Type 


Old type 


Other ID 




(hh:mm:ss.ss) 


(dd:mm:ss.s) 


(days) 


(days) 








V2 


18:24:29.15 


-24:51:06.6 


58 




LPV 


Irr 




V3 


18:24:30.17 


-24:50:18.5 


49 




LPV 


Irr 




V4 


18:24:30.12 


-24:51:36.4 


13.37 


13.462 


W Vir 


W Vir 




V5 


18:24:29.48 


-24:51:52.6 


0.645 


0.644360 


RRab 


RRO 




V7 


18:24:44.19 


-24:50:30.6 




320.0 


LPV, f 


Mira, f 




VIO 


18:24:39.51 


-24:53:26.1 






LPV/Irr 


Irr 




Vll 


18:24:31.48 


-24:51:33.0 


0.542 


0.542767 


RRab 


RRO 




V12 


18:24:43.31 


-24:52:56.1 


0.578 


0.578228 


RRab 


RRO 




V13 


18:24:25.77 


-24:52:33.6 


0.656 


0.654923 


RRab, f? 


RRO 




V17 


18:24:35.84 


-24:53:15.8 


62 


91.7 


RVTau 


RVTau 




V18 


18:24:36.54 


-24:51:49.9 


0.640 


0.640151 


RRab 


RRO 




V19 


18:24:30.84 


-24:51:56.0 


0.3119 


0.335 


RRc 


RRl 




V20 


18:24:33.11 


-24:51:42.9 


0.4948 


0.49774 


RRab? 


RRO 




V21 


18:24:32.94 


-24:51:59.0 




29.93 


LPV 


W Vir 




V22 


18:24:30.98 


-24:52:01.8 


0.323 


0.99538 


RRc 


BL Her 




V23 


18:24:30.25 


-24:52:03.0 


0.292 


0.29231 


RRc 


RRl 




NVl 


18:24:28.89 


-24:52:09.9 


0.748 




RRab 






JN Vz 




-z4:jz:3z.3 


U.z9/ 




KKC 






NV3 


18:24:32.31 


-24:52:32.2 


0.323 




RRc,f 






NV4 


18:24:24.13 


-24:54:36.6 


0.261 




RRc 






NV5 


18:24:27.67 


-24:52:15.0 


0.311 




RRc 




RC133? 


NV6 


18:24:33.45 


-24:52:10.6 


0.308 




RRc 






NV7 


18:24:32.17 


-24:52:00.4 


0.329 




RRc 






NV8 


18:24:33.60 


-24:52:07.1 




13.53 


W Vir? LPV? 






NV9 


18:24:33.05 


-24:51:28.9 






LPV 







5 




Prieto et al.: Variable Stars in M28 (RN), Online Material p 1 
vll 



vl2 





2.0 0.0 



0.5 1.0 

Phase 




vl8 



v20 



- / 

• 


1 


1 

• 


1 






• 


1 


1 

• 


1 


• 

• 

• 


\ 

1 


• 

• 

1 


\ 

1 




GQ 


••• • 


• 

1 


• ••• • 

• • 

1 


• _ 

1 




1 


1 

••• 


1 






• 


1 


1 

• 


1 


• 

t 


1 


• 

1 


1 




> 


.**,••. 

>- • 


• 

1 


.**,••. 

• • 

1 


• _ 

1 




1 


1 


1 






- • • 


1 


1 

• • 


1 


• 

• 

• 
• 


1 


• 

• 

1 


1 






• 

- • 


• 

1 


• 

• •• 
• 

• 

1 


• 

1 



2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

Phase 




• • • • 

\ h K . ft K,. 



0.0 0.5 



_l l_ 



1.0 1.5 

Phase 



Figure 2. Light curves of the ab-type RR Lyrae detected in this work, in the B (top panels), V (middle panels), and / (bottom panels), 
respectively. 
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Figure 3. As in Figure [ij but for the c-type RR Lyrae. 
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Figure 4. As in Figure pi but for the type II Cepheids. 
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Figure 5. As in Figure^ but for the LP Vs. 



